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THE  EXACT  mTH~PLATTMl-MG  PROCEDUilE  IH 
XQNOSPHERIO  PROPAGATIOM  PROBLMS 

ABSTRACT 

Tlxe  exact  earth-flatte-ning  procedtire  prev3.ons]y  developed  fox-  an  iso¬ 
tropic  spherically-stratified  atnosphere^  is  ext:.end3d  to  the  case  of  a 
spherical  earth  and  atmosphere  enveloped  by  a  sharply  boimded  ionosphere.  The 
general  solution  of  the  problem  is  foimulated  as  an  integral  representation, 
f ran  trhich  may  be  derived  eriiser  a  ray-optical  series  or  a  normal  node  series . 
In  the  latter  case,  the  normal  modes  involve  the  normalised  spherical  Kahkel 
function  and  its  derivative.  .An  improved  method  o.C  obtcdning'  the  seres  of 
these  Ihmctions  is  derived  v/Mch  is  riot  of  asymptocic  character , 

.A  spherctdal  geometry  is  .investigated  as  a  basts  for  deali-ng  vdth  prooloms 
of  non-sphericsl  stratificatioa.  Solutions  for  th':  •mgular  on  as  an 

ra'inite  seri  es  of  Be-ssel  funcMons  are  i^ounii,  of  '  ■  >  .cams  tj/pe  r-.=  in  the 
spherical  case  ''Cite  radial  f  i-'’ricn  is  expreoted  .s  a  sua  of  the  im-r.aaliiied 
spherical  Hannet  function  and  ,tt.s  dexv  vett  ve,  the  c.c 'efficients  of  t;.>eae  func¬ 
tions  being  infinite  xerior:  in  teciri.3  -'f  pavers  of  tco  ratio  of  seui -focal 
distance  to  radius,  .'t  is  sbo.m  that  the  aeros  of  I'ue  racial  fun.-j~ion  as  a 
function  of  order',  ri-i-oh  a.i.’e  r-;arired  for  'fc}:.e  .aorme  mode  ::oluti'Oiv  may  be 
foxxad  by  1iie  zbsc.b  p'rc  odU'.e  that  \Ta.s  -ievelcpad  for  +he  spherical  ca'oe 
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VLF  ENHANCEMENTS  AND  HF  FilDEOUTS  DURING 
SUDDEM  IONOSPHERIC  DIvSTURBANCES 

ABSTRACT 

Slnmltaneous  observations  of  short-wave  fade-outs  of  a  3.3o5-*fo/s  signal 
and  sudden  signal  enhancements  of  a  31°15~hc/8  signal  over  substantially  the 
same  transatlantic  path  of  approximately  5400  km  show  no  evident  correlation 
between  the  magnitudes  of  the  two  effects  of  the  SIDo  This  absence  of  correlation 
la  understandable  on  the  basis  of  a  two-layei'  D-regiono 

The  relative  intensifications  of  the  two  D-regioas  will  depend  on  the 
spectral  distribution  of  hard  X-rays  in  the  1-10  A  range  emitted  during  a  flare^ 
which  can  be  expected  to  vary  from  flare  to  flare »  Since  the  increase  in  h-^f 
absorption  is  the  sum  of  the  increases  in  thf.-  tuc  regions;;  while  the  v-l-f 
enhancement  is  occasioned  only  by  thg  changes  at  the  lower  level,  no  correlatiion 
should  result  between  the  two  effects. 

On  the  other  hand,  an  adequate  explanation  of  the  mechanism  of  the  v-l-f 
enhancoDent  is  not  available  on  the  basis  of  present  knowledgeo  Phase  measurements 
show  that  a  definite  decrease  in  height  of  the  lower  boundaijy  of  the  D-regicaa  is 
caused  by  the  flare.  This  reduced  height  causes  reflection  to  take  place  at  a 
level  of  higher  collision  frequent,  which  should  result  in  a  decrease  in  the  effec¬ 
tive  conductivity  of  the  layer  if  the  icaization  gradient  r^aains  the  same.  Conse¬ 
quently,  it  appears  that  an  Increase  in  the  sharpness  of  the  Iowa*  boundary  of 
the  D~region  is  required  during  the  onset  of  a  solar  flare.  The  mechanism  by 
which  this  takes  place  needs  to  be  determined , 
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PART  I 

THE  EXACT  EARTH-ELATTMING  PROCEDURE  IN 
IQHOSFHERIC  PROPAGATION  PROBLEMS 


1,  INTRODUCTION 

In  an  earlier  paper  [l]*,  an  exact  earth-flattening  procedure  was 
given  for  propagation  in  an  inhomogeneous  atmosphere  over  a  spherical  earth 
Thf.s  formulation  led  to  the  realization  of  the  physical  nature  of  the  approxi¬ 
mations  introduced  by  the  usual  earth-flattening  procedure »  In  particular  it 
was  shown  that  the  differential  equation  for  the  height-gain  function  in  the 
usual  earth-flattening  approximation  was  equ.lvalent  to  a  small  change  in  the 
refractive  index  variation  with  height.  In  other  words,  the  physical  problem  is 
changed  somewhat  by  the  earth-flattening  approximation .  The  amount  of  this  change 
or  deviation  Increases  with  height,  but  should  not  be  of  great  consequence  in 
problems  of  tropospheric  propagation- 

In  the  case  of  ionospheric  p?opagation,  the  important  heights  involved  (in 
wavelengths)  may  be  considerably  g^'eater.  Consequently it  appeared  desirable  to 
in’.'estigate  whether  the  exact  earth-flattening  procedure  could  5ji5)rQve  ionospheric 
propagation  analysis.  This  is  on?  ob.iective  of  the  research  conducted  \aider  this 
pa:ft  of  the  contract,  and  is  accomplished  in  See.  2  An  additional  objective  is 
the  extension  of  this  theosy  to  take  into  account  lateral  variations  of  the  re¬ 
fractive  index  (non-horizontal  strati ficeti.oc ) .  For  this  purpose  a  spheroidal 
geometry  is  considered.  This  is  carried  out  is  Sec-  3. 

The  subject  of  ionospheric  propagation,  invol'/ing  complex  layer  distributions, 
isagaeto-ionic  splitting  and  propagation  at  ar-bitraxy  angles  to  the  earth  s  laagnetic 
field,  coupling  between  modes,  etc.,  eacompaBses  aary  ramifications  which  probably 
never  vjill  be  capable  of  a  complete  self-contaiixed  treatment.  Consequently,  for 
•Numbers  in  brackets  refer  to  the  corresponding  numbers  in  the  References  on  p.  29. 
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purposes  of  the  present  stud7  we  shall  adopt  an  often-used  idealization  of  the 
ionosphere  in  order  to  confine  attention  to  the  specific  objectives  stated  above. 

Sbr  this  purpose  the  ionosphere  will  be  considered  to  be  sharply  bounded  and  of 
unlfonn  electrical  properties.  This  assumption  is  the  one  usually  made  in  study¬ 
ing  v-l-f  ionospheric  propagation,  so  that  the  results  will  be  of  chief  interest 
in  this  frequency  range.  It  la  then  logical  to  consider  only  a  vertical  dipole 
source,  since  this  is  the  only  effective  form  of  radiator  at  these  frequencies. 

2.  SPHffilGALLY-STRATIFIED  IONOSPHERE 

A  rigorous  formulation  of  the  field  due  to  a  vertical  electric  or 
magnetic  dipole  in  an  inhomogeneous  isotropic  atmosphere  over  a  spherical  earth 
was  given  by  fM.edman  [2].  For  plane  geometry,  this  was  extended  by  Walt  [3]  to 
Include  the  essential  mixed  polarization  effects  due  to  the  anisotropy  of  a  sharply 
bounded  ionosphere.  For  completeness,  a  rigorous  formulation  of  the  spherical 
problem  (with  a  ahu.'rp  ionosphere  boundary)  will  be  sketched  here,  Thi.s  formulation 
will  be  given  in  a  form  adapted  to  direct  introduction  of  the  earth-flattening 
procedure. 

In  the  isotropic  case  treated  hy  Friedman,  it  is  possible  to  formulate  separate¬ 
ly  the  cases  of  vertical  electric  and  vertical  magnetic  dipole  sources,  corresponding 
to  vertically  and  horizontally  polarized  fields,  respectively.  In  each  case,  the 
various  field  con^wnents  are  derivable  from  a  Hertz  vector  whose  direction  is 
radial.  Actiially  this  Hertz  vector  (within  an  appropriate  multiplying  factor)  is 
nothing  more  than  the  radial  component  of  the  electric  (magnetic)  field  in  the  case 
of  the  radial  electric  (magnetic)  dipole  source,  since  all  other  components  are 
derivable  from  the  radial  coii5X)nent8  (see,  for  example,  Schelkunoff  [4])o  In  the 
anisotropic  case,  however,  electric  and  magnetic  modes  are  coupled  in  the  ionosphere, 
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so  that  the  problem  must  be  formulated  in  terms  of  mixed  components  from  the 
outset, 

2ol  Formulation  of  the  Problem 

The  geometry  of  the  problem  is  shown  in  Fig,  1,  A  vertical  dipole  of 
(infinitesimal)  length  i  and  current  I  is  located  at  R  =  b,  the  boundaries  of 


Ionosphere 


Vertical 


Figc  1  -  Geometry  of  spherical  earth,  with  concentric  sharply- 
bounded  sfiherical  ionosphere,  excited  by  dipole  source, 

earth  and  ionosphere  being  at  R  =  a  end  R  =  h,  respectively- 

Consideration  of  the  physics  of  the  problem  will  assist  a  proper  formulation , 
Thus,  the  primary  field  due  to  the  source  will  give  rise  to  a  field  which  has  a 
polarization  determined  by  the  direction  of  the  source  current.  This  primary  field, 
in  turn,  will  give  rise  to  reflected  components  at  the  boundaries  of  the  earth  and 
ionosphere.  The  ionosphere  will  introduce  magneto-ionic  splitting,  so  that  new 
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polarization  components  will  arise  there.  Ptom  these  facts  it  is  cleajf  that  a 
combination  of  electric  and  magnetic  Hertz  vectors  must  be  used  for  derivation 
of  the  fields.  The  two  conqponentS)  in  generalf  will  differ  in  amplitude  and 
phase,  SO  that  we  must  represent  the  radial  Hertz  vector  by  a  column  matrix  of 
the  foxm 

Cs]  = 

where  is  the  unit  radial  vector,  and  the  subscripts  e  and  m  refer  to  electric 
and  magnetic  modes,  respectively. 

Consider  first  the  electric  component  He  and  write  it  as  the  sum  of  a 


primary  and  a  secondary  field 

IIq  =  Hi  +  IIj. 

(2.2) 

Now  put 

n, 

(2.3) 

Hi  is  stimulated  by  the  vertical  source  current,  while  IIj  arises  firom  reflection 
at  the  boundaries.  Ihen  the  corresponding  fields  are  derivable  from  the  equations 


E<  *  kR(p,tP,>+^  3Md[^{R(P,  +  P3i}],  (2.4) 

providing  that  F^,  is  a  solution  of  the  Inhomogeneous  reduced  wave  equation 

V*P,  +  =^0-  (2.6) 

and  F^  is  a  solution  of  the  homogeneoxis  equation 

V'Pjtk'Pj.O.  (2.7) 
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First  consider  (2=6) «  The  current  density  aaay  be  related  to  the  dipole 
jnoiBent  II  by  Integrating  over  the  source  regions 


<p  2  ft 'll 

11  [  j  JiR'^sin$  dddfdf?^ 


so  that 


C  d  o 


li 


EttR*sin^ 


S(S)SCR~b). 


=8) 


Sin*e  the  right-hand  side  of  (2o6) ,  in  virtue  of  (2=8),  is  zero  every</here 
outside  the  point  (bj0,0),  the  solutions  of  (2  =  6)  can.  be  assembled  froia  solutions 
of  the  corresponding  homogeneous  equation 

+  k^P,  =  0.  (2=9) 

Hence  we  can  separate  Pi  in  the  fora 

R  P,  •-  T(&)  L),  (R)V(<p)  (2=10) 

where  T,  U,  and  V  are  functions  only  of  8  R,  ard  .p,  respectively =  (2^9)  then 

separates  into  the  equations 


de^ 


cot€ 


\ 


T 


O, 


(2=11) 


dR* 


ilY. 

dcp^ 


+ 


V 


(2=12) 


^2=13) 


in  ’.rhlch  s  and  it  are  the  separation  constitots,  wbj.cb  as  yet  are  arbitrery,  and 
ultims.tsly  will  be  fixed  by  the  boundary  ■onditj.or.s ,  The  various  solutions  of 
(2-9}  are  characterised  by  different  '-^alues  of  s  snd  m,  including,  possibly,  complex 
values  = 

ds  tab.8  a  to  be  an  ^Xiteger  '"o  ovnor  •■•at  V  htVE;  .Sn-periodicity  in  and  write 
the  Ro3.uticn  of  (2=13)  in  the  form 

\i^  ~  cos(r^l^■^ 


(2=14) 


Consequently  all  solutions  of  (2o9)  with  Sit-periodioity  in  <jp  may  be  obtained 
from  the  representation 

RP|*Z  fA(6)T(0)U,(R)V^rq>>^s,  (2.15) 

m«Oj( 

where  the  amplitude  function  A(8)  and  the  path  C  ia  the  complex  s-plane  are  as  yet 

unspecified.  In  general,  C  will  extend  over  an  infinite  range.  A(s)  and  C  rnty  be 

det«nnined  integrating  (2.6)  around  an  infinitesimal  region  enclosing  the  dipole 

sotirce.  It  can  be  shown  that  A(s)  =  s  and  /  =  /  »  provided  that  T(0)  =  1,  so  that 

RP,  =£  /sTU,V„d5,  (2-16) 

m«o  o 

where  T  is  a  solution  of  (2.11),  Ui  is  a  solution  of  (2.12),  and  is  given 

(2.U). 

2.2  THa  Angular  Function  T 

In  [1]  it  was  shown  that  a  solution  of  (2.11)  for  m  =  0  is 

T=  i  (2-17) 

n«0  * 

where 

in  which  2^(80)  ia  a  cylinder  function.  In  order  that  (2.17)  have  the  property 
T(0)  =  1  as  required,  we  must  choose  the  cylinder  function  to  be  the  Bessel  functicni 
J^(s0),  and  =  1.  Consequently  the  required  solution  of  (2.11)  for  m  =  0  may  be 
written  as 

Z  tt„vse>"  J„(»0).  (2.18) 

n*0 

It  may  be  shown  that  a  lower  bound  for  the  absolute  convergence  of  (2.18)  is  j9|  =  2, 
so  that  this  covers  a  sector  greater  than  +  x/2. 

We  now  extend  this  type  of  solution  1x>  the  case  m  0. 

Introducing  the  new  independent  varifible 

X  =  se,  (2.19) 
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and  denoting  the  dependent  variable  by  y,  (2 oil)  becomes 
<**4- ictsy'  +  (l  -  =  O. 

We  write  (2.20)  in  the  form 


(2.20) 


where 


O-)!  »  2*^  B|ib/(2^)! « 


(2.21) 


(2.22) 


the  being  the  Bemouilli  numbers. 

Assuming  a  solution  of  (2.21)  of  the  form 


(2.23) 


we  obtain 


L(«^  * Z  s-^LiC*^)  *  Z  Z  .4- (2|i-i)w»x-*^J  • 

»  V>0  ®  w«0  fist 

equating  coefficients  of  like  powers  of  8>  we  obtain  the  systofi  of  equations 

L(i^)*0, 


(2.24) 


0  a  • 


A  solution  of  the  first  eqiiation  is 

7o  =  Ztt(x), 

where  Zj^  is  any  cylinder  function.  The  second  equation  then  becomes 

L(ya)  =  ai(xZ^+m*Zjn) 

=  ai[n(m<'l)Zia-xZjj<.i] 

hy  introducing  the  function 


*  A^Zin+n  ^X)j 


which  has  the  property 


(2.25) 


(2.26) 


(2.27) 
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(2  .,25)  becomes 


L(ya)  =  ai[m(mH)Cm,o-Cni,i3o 

Now  using  the  property  (2.27),  the  solution  of  this  equation  is  seen  to  be 

£|y  induction,  we  infer  that 

V 

Hence  the  solution  of  (2.21)  should  be  expressible  in  the  form 

•  n»0  A«0 

The  following  recursion  formtalas  for  *^m,n  aa’e  easily  obtained  fi^am  the 
recursion  formulas  for  the  cylinder  functions: 


(2.28) 


(2o29) 


(2.30) 


~  *"A,P  ^m,n*A+p  » 


where 


-  /  ■/  Z*~ 

pi  (^-p)!(m+n+p)j  ' 


(2.31) 

(2.32) 


(2.33) 


If  we  substitute  (2.30)  into  (2,23),  and  use  (2.31)  land  (2.32)  to  eliminate 
powers  of  x  on  the  right-hand  side,  we  obtain 

^2(n+i)A„,,,C„^rt  -  ^  ^WiP+H+J  ^»i,p+p+j+i)» 


where 


Ao  =  l. 

D  (2p+l)m*+m+Zp. 

Ely  equating  coefficients  of  like  orders  of  the  function  Cin,n  on  the  two  sides  of 


this  equation,  we  obtain  tb.e  recursion  formula 
I  fA  *<<• 

An+I~  9ly,J.n  |2r  2-,,  0  Cpu  ^-p-B 


=  DC  -7  7  c  U  u. 


(2,34) 


Consequently,  the  required  solution  of  (2.11)  is 
T*  i  A«(50)"CL.„(58). 


{2o35) 


8 


The  advantage  of  using  an  expansion  for  T  in  texiss  of  Bessel  functions, 
instead  of  the  standard  expression  in  ternis  of  the  associated  Legendre  functions, 
is  that  a  nore  accurate  calculation  is  possible  than  by  the  use  of  the  asymptotic 
expansion  for  the  latter  functions. 

2,3  The  Radial  Puncticm  P 

With  T  as  given  by  (2.35)  the  solution  of  (2.9)  is 

I  A  „ (ssy  (sS) cosiin^*ym)  U  (2,36) 

where  Aq  ^  1  and  Aq  is  given  by  the  (2,34). 

The  integral  along  the  positive  real  s-axis  in  (2,36)  may  be  transformed  into 
an  Integral  along  the  entire  real  axis  In  the  following  way: 

Write 

and  note  from  (2,11)  and  (2,12)  that  T  and  U  are  even  functions  of  s.  In  the 

integral  corresponding  to  make  the  substitution  s'=se"’*’'^  ,  whereupon 

the  integral  for  that  term  becomes 

o 

in  view  of  the  fact  that  the  integrand  is  an  even  function  of  s.  Then  (2,36) 
becomes 

R p,  =  i  Z  2  r A,  (sfi)"  (se)  cos  Cw<p+ -4,,)  U,  s  ds,  (2,37) 

This  form  is  adaptable  to  evaluation  by  residues  or  by  stationary  phase,  depending 
on  whether  a  normal  mode  representation,  or  a  representation  in  terms  of  rays  is 
desired. 

The  function  Uj,  is  to  be  fixed  by  the  boundary  conditions.  These  require 
that  the  tangential  electric  and  magnetic  fields  be  continuous  at  R  =  a  and  R  =  h. 
For  this  purpose  both  the  electric  and  magnetic  components  of  [g]  will  be  required. 


(2,38) 


(2.39) 

(2.40) 

(2.41) 


Hence  we  now  consider  the  magnetic  component  Hq  in  (2.1) »  and  write 

a.  = 

Then  Pg  satisfies  the  homogeneous  equation 

V*Pe  +  k»Pj*0. 

The  corresponding  fields  then  are  derivable  trem  the  equations 

e„-curl(gP.), 

Solutions  of  (2.9)  and  (2.39)  may  be  written  in  a  form  similar  to  (2.36)  as 
follows) 

RPa^SS  f  (2.42) 

iniOii«o  ^ 

RP3  *  £  i  f  5  An  (s8)  co5(m<f 't-r'm)U3  s  ds.  (2.43) 

n*o  JJj 

The  constants  S  and  5  are  to  be  determined  by  the  boundary  conditions  at  R  »  a 
and  R  =  h. 

Corresponding  to  the  physical  picture  of  reflection  at  the  boundaries »  ve 
expect  a  mixture  of  upgoing  and  downgoing  waves  in  the  region  a<R<h.  We  then 
pick  the  two  independent  solutions  of  (2.12)  to  correspond  to  upgoing  and  downgoing 
wavesy  and  denote  these  by  and  f  respectively.  A  similar  choice  is  made 

for  Ug  and  U3.  The  total  field  in  the  various  regions  then  can  be  derived  from  a 
radial  P  function  which  has  the  matrix  form 


s[p]  -  5  [’’'4”]  • 


(2.44) 


in  which 

RPq  =  ^n^^n»  ®  “  l»2f3e 

The  boundary  conditions,  being  Independent  of  6  and  9,  lead  directly  to  the 
statements 

=  Tg  =  Tj, 


Vi  =  Vg  =  Vg. 
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Nov  ve  put 


UjL  =  ^ 

=  B.<‘>  +  V®'i 


Il3  =  BjW 


+  Bj<»>.  J 


and  introduce  the  reflection  coefficient  at  the  ground 


(2.45) 


(2.46) 


where  f,  and  are  the  reflection  coefficients  for  vertical  and  horizontal 
polarization,  reapectively.  Then 


urta)  «  fiUj'ctt), 

U8*’(d)sp,uf£a). 


(2.47) 


At  the  ionoeidiere  the  reflection  coefficient  is  a  tensor 

so  that 

Uf(h)]  +  p,aUf(h), 

Ui'Yh)*  pt,[U,'’’(h)'t  U3**(W]  + 

Jlnally,  at  R  =  b  we  have  the  discontinuity  condition  for  the  first  derivative 
of  Ui  in  terns  of  the  dipole  mcnient  [2] 

^1"-'’“  .  _ 
while  U],  Itself  is  continuous  at  R  =  b. 

The  radial  functions  U2,  Uj  satisfy  the  same  type  of  differential  equation  as 
U],,  i.e.  (2.12)0  If  we  denote  the  two  Independent  solutions  of  this  equation  by 
and  u^*^,  respectively,  where  u^^^  represents  a  downgoing  wave  and  u^®^  an 
upgoing  wave,  then  we  may  write  in  the  various  height  regions 


c 

II 

b<R< h , 

(2o49a) 

U,  =  lSau‘« 

fl.<R<  b  t 

(2o49b) 

U2  * 

o<R<h, 

(2o49c) 
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+  If,  <KR<h .  (2.49CI) 

The  boundary  conditions  then  yield 

d&a/ei  >  pt  u‘’*(a)/u‘*V),  (2.50a) 

9t/€t  *  Pi  u‘'‘<ft)/u**‘(a>,  (2.50b) 

=  f’l  u“‘to)/u‘*»(fl),  (2.50c) 

St  u^'lh)  *  [  |.«  (e^  t  4  Jbi]  u‘"Uh)  (2. 50d) 

Uiu‘'’(h)-  [p,,(».^g^  +  (2.50e) 

»,u‘" '(fa)  s  tf,  uV»(b)  4  u«»'(t>)  +  K,  (2.50f) 

=  ^,u“*(b)  f  (2.50g) 

The  seven  equations  (2.50a-g)  are  sufficient  to  detennine  the  seven  constants 
®*f  "Si*  'Cft  -6,*  Sg,i  »  tt'-  They  are  given  by 

•  K/[(i^b-  u5*>'(b)]  *ilKu5«(b),  (2.51a) 

(2.51^) 

®«*  (p'Va+ W*'*  (2.51c) 

(2.51d) 


where 


Pti  - 

f*'^K 


7b  ’ 

■  u‘'’(h) 

u‘*)'(b)  ’ 

_ _ _ 

CpM  -  (I,  -  fu  f.)  <^8 

A  =  p„  p»j,  —  p^, 


(2o51e) 

(2<>51f) 

(2.51g) 

(2.52a) 

(2.52b) 

(2.52c) 

(2.52d) 
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primes  denoting  R-derivatives  evaluated  at  the  argument. 

We  XKJW  evaluate  the  form  of  the  radial  functions  u^^^  and  u^^^.  These 
are  solutions  of 


u"+  (k*-'-|i)u«0.  (2.53) 

The  solutions  of  this  equation  corresponding  to  downward  and  upward  waves  are 
the  nonnalized  spherical  Hankel  functions  [5] 

u‘'>  .  h‘,'’(kiy  •  “(kW,  (2.54) 

u®  =  h'"(l<IO  -  (2.55) 


respectively,  where 


(2.36) 


With  those  functions  inserted  in  (2.49),  the  expressions  (2.37),  (2.42),  (2.43) 
give  the  values  for  RP^  in  the  space  a^<h,  frca  which  the  fields  may  be  evaluated 
Py  (2.4),  (2.5),  (2.40),  and  (2,41). 


2.4  Evaluation  of  the  Integral  Representation 

Two  different  methods  are  available  for  evaluating  the  integral  expres¬ 
sions  for  RP.  the  method  of  stationary  phase,  the  result  may  be  expressed  as 
a  sum  of  rays  reflected  alternately  a  number  of  times  from  the  ionosphere  and  the 
ground.  Sy  the  method  of  residues,  on  the  other  hand,  the  result  is  obtained  as 
a  sum  of  normal  modes,  or  waveguide-type  waves.  We  shall  investigate  the  latter 
type  of  solution  in  order  to  bring  out  the  fact  that  the  approximations  usually 
made  actually  change  the  physical  problem  from  that  of  a  homogeneous  atmosphere 
to  that  of  a  slightly  inhomogeneous  atmosphere. 

Since  the  coefficients  in  the  integrand  (®3  “  involve  the  y-functions 
defined  above,  which  are  ratios  that  are  functions  of  s,  the  integrand  has  poles 
at  zeros  of  th^  denominator  in  these  ratios.  Consequently,  if  we  deform  the 
integrand  from  the  original  contour  along  the  real  s-axis  into  the  appropriate 
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half  of  the  complex  plane,  lixe  integral  may  be  evaluated  in  terms  of  the 
singularities  of  the  integrand  in  that  half-plane c  In  addition  to  the  poles 
j\ist  mentioned,  there  is  also  a  branch  point  vrhere  the  order  of  the  spherical 
Hankel  functions,  p,  is  zero.  This  can  be  seem  from  (2.56).  This  has  branch 
points  at 

®  ~  i  ^i, 

The  integrand  vanishes  at  Infinite  values  of  s  in  the  lower  half -plane. 
Coxisequently  the  integration  path  is  deformed  into  the  contour  shown  in  Hg.  2. 

The  Integral  then  is  the  negative  sum  of  the  residues  at  the  poles  in  the  lower 
half-plane,  plus  an  integral  aroimd  a  branch  cut  along  the  negative  imaginary 
axis  from  -l/2.  Friedman  [2]  has  discussed  the  importance  of  the  branch-cut 
Integral  and  has  shown  that  it  is  negligible  in  practical  cases.  Wait  [3],  on 
the  other  hand,  attempts  to  avoid  the  branch-cut  integral  by  making  a  double  tra¬ 
verse  in  the  lower  half-plane,  but  his  procedure,  in  effect,  is  equivalent  to 
neglecting  this  integral.  This  integral  represents  the  effect  of  the  currents  which 
penetrate  into  the  ground,  and  thus  is  essentially  a  part  of  the  ground-wave  field. 
In  the  case  of  a  perfectly-conducting  ground  the  integral  vanishes  altogether. 

The  matrix  ^P]  in  (2.44)  has  an  integral  representation  vdiich  can  be  assem¬ 
bled  frem  (2.37),  (2.42)  and  (2.43)  by  using  the  U-functions  given  in  (2.49).  Poles 
of  the  integrand  are  those  of  the  functions  and  M.  The  principal  poles  of  in¬ 
terest  in  detennining  the  normal  modes  are  those  of  M.  The  investigation  of  these 
poles  is  a  separate  probloa  in  its  otoi  right  v/hich  we  shall  not  go  into  here.  The 
poles  of  Pa,  since  pa,  \iltimately  can  be  e:3q)ressed  in  teuas  of  y-functions  and  the 
properties  of  the  reflecting  medium,  can  be  expressed  in  terns  of  the  two  limiting 
cases  pi, and  pi,  similar  to  the  way  in  which  Bremmer  [6]  treated  the  tropo- 

s|heric  case .  These  can  be  detenoined  from  tb.e  ceros  of  and  Cb),  respec¬ 

tively.  Thus  we  consider  the  method  used  for  the  determination  of  these  seros. 
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2o5  The  Complex  Zeros  of 

The  zeros  of  and  are  the  same  as  those  of  Hp^®^  and  , 

These  are  foimd  by  the  Debye  method  of  steepest  descent,  and  are  usually  ex¬ 
pressed  in  tents  of  Airy  functions,  or  Hankel  functions  of  order  one-third  o  The 


procedure  is  to  vnrite 


h“‘<z)  »  ^  P(w-Tr/a)3 dwif  dw, 


(2o57) 


w*  "Z 

expand  the  exponent  F(v)  in  a  Taylor^s  series  about  the  point  where  F"(v)  =  0, 
and  draw  the  contour  Vq  so  as  to  pass  through  the  two  points  (stationary  points) 
at  which  F'(w)  =0.  Dy  truncating  the  Taylor'^s  series  expansion  of  F(w)  at  the 
third  derivative  term,  we  obtain 


Since 


we  have 


F(w)js  F(w*)  +  (w-w.^F'(w.)  + 


F"(We)  *-iaeos\«4,  »  O, 


Wq  —  x/2, 


F(wo)  =  0, 

F'(wo)  =  -i(z-p), 


F«'(w^)  =  iz. 


Consequently,  upon  putting  v-Vq  =  u,  (2 <,57)  becomes 
Hp  (z)«>5pje  du, 

where  the  contour  Ug  is  merely  shifted  to  the  right  x/Zc  A  sisqple  change 


of  variables 


t=  (iz/zy^u 


(2o58) 
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results  in 

(2.59) 

where  the  contour  Iig  in  the  t-pl@ne  is  shown  in  Pig.  3.  The  integral  in  (2.59) 
may  be  expressed  in  terms  of  the  Airy  Functions,  or  Modified  Hankel  Functions  of 
order  one-third  [7].  Using  the  notation  for  the  latter, 

u* 

we  obtain 


^p\z) «  hz  («;). 


(2.60) 


t-  plane 


Pig.  3  —  Contour  for  Modified  Hanxel  Functions 


Then  from  (2,55) 

Consequently  the  zeros  of  63(5)  (tabulated  in  [7])  give,  in  first  approximation, 

the  zeros  of  and  of  hp^^Hz). 

It  was  pointed  out  in  [1]  that  the  approximation  (2.61)  is  equivalent  to  a 
change  in  the  physical  problem.  This  is  immediately  evident  from  the  fact  that 
hp^2>{z)  la  a  solution  of  (2.12),  while  h2(i;)  is  a  solution  of  Stokes*  equation 

4%  t  O.  (2c62) 
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The  physical  problem  corresponding  to  (2o6l)  may  be  found  as  follows: 
We  first  &«;parate  (2,9)  by  writing  it'i  the  fom 

P|  =  T(e)Uo(R)V(«p>, 
wherety  the  radial  equation  becomes 

dR^  R  dR  R»/ 

How  by  introducing  the  transformation 

=  a  lo^(R/o.), 

the  radial  equation  becomes 


d»a  1  du.  ^ 

^  dr\ 

Next,  putting 

Uo  » 

we  obtain 

dy]^ 

where 

=*Ja  (**•'■  ’-4)  = 

iw) 


(2,63) 


To  reduce  this  to  Stokes'  equation,  we  must  have 

k^e^**/*  *  k*(i  +  <5i?),  (2,64) 

where  Icq  and  q  are  constants.  It  is  evident  from  (2,63)  that  q  -  2/a  in  order  to 
satisfy  the  equation  for  small  In  this  case,  if  we  put 

s*  ^  o-yki, 

^  =(ko/<^^(f-S*+c^^U 


we  obtain,  finally,  Stokes’  equation 

el*u.a  _ 

d  ■5*-  ^  o, 

a  solution  of  which  is 


He  = 

In  order  to  arrive  at  this  solution,  however,  k  must  be  a  f met  ion  of  which 


IS 


satisfies  (2.64)  1  In  tema  of  the  variable  R,  this  requires  that  k  have  the  form 

(2.65) 

If  we  put 

R»  a-f  H  *  a(l+ 


then  ve  have 


yt(2) 


i7? - 


4* 


Thus  the  refractive  index  (2.65)*  corresponding  to  Stokes'  equation,  decreases 


aonotonlcally  with  increasing  height  H  above  the  ground  level  a,  whereas  the 


original  probl»i  dealt  with  a  constant  refractive  index. 

To  obtain  a  higher  order  approximation  for  the  zeros  of  one  Jia7 

follow  a  procedure  due  to  Olver  [8]  and  Chester,  Friedman,  and  lirsell  [9],  whereby 
a  change  of  variable  is  introduced  so  that  F(w)  in  (2.57)  becomes  precisely  the 
exponent  in  the  Integral  of  (2.59) t 

F(w)  *  t’/j, 

dw  C»t*  ^  9  + 

dt  *  F'M  -a(»*I»iw-p)  ’ 

Then  (2.57)  becomes 

i  / e'«*^ dt.  (2.66) 

Au 

Ejjr  expanding  ^  in  a  double  series  of  the  form 

^  -  ip-cftsr  +  f  c)„tr(f*«r,  (2.67) 

m»0  m»e  ' 

and  Integrating  (2.66)  termwise,  an  asyiiq>totic  esqjanslon  is  obtained  in  terms 
of  ha(?)  and  l4(^) 

h‘*\z)-  8„>,  (2.68) 

‘  m»* 

in  which  the  coefficients  Ag^  and  Bq  involve,  in  general,  inverse  fractional  powers 
of  z,  except  that  Aq  =  1. 

The  above  pzncedicre  is  a^n^totic  because  the  series  expansion  (2.67)  has  a 
radius  of  ccmver^ence  which  is  limited  h7  the  next  zero  of  f'(w),  which  occurs  at 
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w  “  it-sln"'{p/z),  while  the  intenreil  of  integration  extends  to  infinity* 

An  alternative  evaluation  of  (2*57)  which  is  not  of  asymptotic  character 
nay  be  developed »  however*  This  does  not  appear  to  liave  been  reported  previously* 
We  write 

F(w)  -  ¥^/s  +  Si, 

where  ^  and  t  are  given  by  (2*53),  and 

»»  £  t”, 

nM 

i\*C  t\*S^ 

Next  we  write 


+  ^TSi 


and  expand  e^  in  the  absolutely  convergent  series 
m*o 

*p(t)«  i  +  i 

tii«S 

The  integral  (2.57)  then  beccanes 

Ua 

Termwise  integration  then  yields  an  expression  of  Identically  the  same  form  as 

(2*68), 

Hp^Cx)  ‘ 


(2*69) 


where 


=<(»;)  =  1  + 

m^S 

Ajg(5)  and  %(^)  being  polynomieLls  in  j;  of  degree  m/2  or  less* 

Then 

C<x)  = 

=  (a4)''*i>*z'*e"‘’^'”*{<<('4))ij(?)  +  /3ft)hatoi.  (2.70) 

It  is  evident  from  (2,69)  that  the  zeros  of  H^^(z),  for  given  z,  differ 
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slightly  from  the  zeros  of  hg(^),  where  5  is  related  to  z  and  p  by  (2»58)o  In 
order  to  find  the  values  of  p  for  which  (z)  is  zero,  we  can  proceed  as  follows: 
Denote  the  zeros  of  hjg(^}  hy  so  that 

ha(5o)  =  0.  (2.71) 

From  (2.69)  we  then  find 

H“\z)  » 

The  value  of  p  correspojoding  to  ?o  ^  ^p^^  (’b)*  denote  this  zero 

by  (T*!  1. o Q« } 

e^Hz)  =  0, 


and  put 


p  =  <r  -  q.  (2.72) 

(Tt  which  as  yet  is  unknown,  corresponds  to  a  value  of  4  which  we  denote  hy 

'S  *  'i)  -b  ’?i«  (2.73) 

so  that  is  small  compared  to  Then  (2.69) 

«  0  »  0(('5o+  «i)  t  /3(Se.+>5i)  h'z 

Ve  now  expand  of,  /3,  hzf  ^  Taylor's  series  about  :;o,  and  make  use  of  (2.71): 

«<o+  **’  •» 

.s(«5)«  |5o+  'S.|3o+ 

3  _ 4  y 

hz('S)- 

*  fi* 

h4^0*  +'3o^-^+- '  •} 

Consequently  we  obtain 


>7,  +  ;7a/S(i5)  ~0. 

This  is  a  series  in  'ft,  whose  coefficients  are  known.  The  value  of  f j  then  may  be 
obtained  ty  successive  approximations.  Then  from  (2.58),  (2.72)  and  (2.73), 
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we  obtain 


<1  s  .  (2,74) 

The  reqiiired  zeros  of  (z)  then  are 

p  =  <r-  q. 

These  likewise  are  the  zeros  of  the  loodified  spherical  Hankel  function  h^^{z) 
given  bgr  (2.70). 

The  detailed  results  obtainable  this  procedure  will  be  reserved  for  a 
later  investigation, 

3.  WOM-SPHERICAIJ.T  STRATTFIED  I0H03PHKRE 

In  the  treataent  in  See.  2,  the  earth-ionosidiere  region  was  assumed  to  be 
spherically  syanetrical  (so-called  "horizontally-stratified”  medium).  This 
situation  is  not  strictly  true,  in  general ,  so  that  the  above  type  of  analysis  is 
an  idealization  which  should  be  considered  as  only  a  first  approximation  to  the 
tzxie  state  of  affairs.  For  example,  there  are  situations  of  practical  interest 
where  the  reflecting  layers  are  tilted  with  respect  to  the  horizontal . 

In  order  to  Introduce  a  form  of  non-sphericel  stratification  which  may  be 
applicable  to  such  situations,  we  consider  the  case  of  a  spheroidal  gecsnetx^’', 
where  the  earth  and  ionosphere  6tre  coordinate  surfaces  of  a  family  of  spheroids, 
either  oblate  or  prolate  in  fora.  We  give  below  the  extension  of  the  exact  earth- 
xlattening  procedure  to  this  non-spherlcal  gtxsnetry, 

3,1  Formulation  of  the  Problan 
The  reduced  wave  equation 

V^P  +  k*P=  O 

may  be  separated  in  spheroidal  coordinates  into  radial  and  angtilar  differential 
equations  as  in  the  spherical  case.  For  the  oblate  spheroid,  these  are 
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^)U  =  0, 

^  +  m^V  =  O, 


(3ol) 


where  f  is  the  semi-focal  distance*  and  a  typical  space  point  has  the  rectangular 
coordinates 

X  *  f  coshl^cosdcos^p, 
y  =  f  cosh  f  cos6  s  in  9, 
zs  fslnh^sin®. 

The  corresponding  equations  of  the  prolate  spheroid  are 


-  ;^)U  -  O, 

cote  ^  ..  (kf  Sm'fi  +  =  O.  (3.2) 

m*V=0, 

for  which  a  typical  space  point  has  the  rectangular  coordinates 

x*  f  sinh^sififi  j 

y-fs'iMhl  sm0  sincp, 
a  *  f  cash  §  c<»»  6, 

We  shall  treat  the  oblate  case  in  detail*  since  a  comparison  of  the  correspond¬ 
ing  eqviations  of  (3<.l)  and  (3o2)  shows  that  a  change  from  (3«>1)  to  (3»2)  can  be 
effected  by  simple  transformations, 

3c2  The  Angular  Function  T 

We  consider  first  the  angular  function  T,  Introducing  the  new  indepen¬ 
dent  variable  x  =  s©,  as  in  Sec,  2,2,  the  second  equation  of  (3ol)  becomes 

^  *5  “‘S)  ^ c«‘{D]t=  O.  (3.3) 

This  may  be  written  as 


(3.3) 


L(n3T"+|T'+  [■jr-icotg)]T'-[m‘[Ji-.^csc*^)]  +  )^.m'g)}T 

=  I,  +  iVC-Tljl)*"  (3.4) 
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whcrtii 


tc  «  kf, 


(3.5) 

(3.6) 


z*^6jLA2}0!  (3.6) 

(3.7) 

(3^4)  is  wiiiin«T»  In  fcnm  iijo  (2.21),  and  dlffsrs  fraii  it  (sily  in  the  presence 
on  l^e  right-hand  side  of  the  itdditianal  term  b,i.K:^(i</s)*^  T  o  This  term  has  the 
S8isi>  power  of  z  as  the  (p+l)-t|;irm  immediately  preceding  it  in  (3.4).  Consequently 
we  isan  imaediately  write  the  e(|;j!ilution  of  (3.4)  as 

T.^^0:fOm,n,  (3-8) 

the  eoeffielent.  (1„  ue  |  iTeD  by  the  recursion  foniula 

,  tt-l 


2(«+!) 


Theii'efore  the  form  of  solution  g|i^sQ  iJ“  (2e35)  is  directly  applicable  to  (3.3) > 

I 

which  thus  has  the  solution  I 


T*  X  (3.10) 

ni*o  I 

3o3  The  Radial  Function  U  | 

We  now  consider  the  rMial  function  U,  which  .satisfies  the  first 

equation  of  (3.1)  o  Our  aim  will  ^  obtain  a  solution  of  this  equation  similar 

■( 

to  that  found  in  Sec.  2,3.  Then  fields  will  be  obtaiiable  in  terms  of  an 
integral  representation  of  the  ftJk  give°  in  (2.36).  We  ishall  be  interested  in 
the  aorasal  mode  solution,  which  is'  obtainable  fpcam  the  reijidues  of  the  integral 

representation.  These  residues,  in  the  spherical  case,  ultimately  may  be  based 

,1 

on  the  2WOS  of  the  function  U  wMileh  represents  an  upgoini;?  wave,  and  its  first 
derivative.  Consequojtly,  we  shal-H  seek  solutions  of  the  radial  equation  similar 
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to  those  given  in  (2o54)  and  (2.55)  for  the  spherical  ease,  and  then  will  investi¬ 
gate  their  complex  zeros  as  a  function  of  order. 

The  radial  equation  in  question  is 

^  t  -a  o.u) 

We  sedc  to  cast  this  into  a  fozm  which  resembles  the  spherical  equation. 

We  first  note  that  the  transition  to  the  spherical  case  is  effected  by  alloving 
f  coshC-^R  as  f->Q.  Hence  we  are  .ed  to  Introduce  the  change  of  independoit 
variable 


and  the  new  dependent  variable 

Then  (3*11)  becxmes 


f  coshC  ~  Rf 


u  s  RU. 


al*u 


■f  — +_l_  rk.R*  -5*+ =  o. 

R(R*ufa)dR  J 


Next|  we  put 


z  -  kRy 
a  =  kf , 


whereupon  (3.13)  becomes 


»(a*-a*)  ^**-a*V‘  z*  r 

primes  denoting  derivatives  with  respect  to  z. 

To  eliminate  the  first-derivative  tezm,  we  put 


m 


V; 


(3*16)  then  is  replaced  by 

r.S*-A*  yw  c  o 

We  now  rearrange  (3.17)  in  the  form 

(I 


(3.12) 

(3.13) 

(3.U) 

(3.15) 

(3.16) 

(3.17) 
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I 


(3ol8) 


where 


, .  m*-  Hi 


Jt-|cti-o],  n£a. 


(3ol9) 


-[i+  1^2:  a. 

If  the  right-band  aide  of  (3<»18)  were  zero  (i.e.,  a  =  0),  solutions  of  the 
equation  would  be  the  nozinalized  spherical  Hankel  functions  gives  by  (2o54)  and 
(2.55) «  Hence  we  seek  a  sisiilar  solution  to  (3.18).  Solutions  of  the  radial 
equation  as  a  series  of  solutions  of  the  spherical  Bessel  differential  equation 
are  available  in  the  literature  (see,  for  eocample,  [10]),  but  we  shall  find  it 
more  ccmvenient  to  deduce  directly  a  special  form  which  is  suitable  for  the  nozmsO. 
mode  problemo 

The  form  of  (3*18)  suggests  a  series  solution  in  a/z.  Such  a  solution  may  be 


formed  in  the  form 


•>>)].  (3-20) 

where  is  a  solution  of  the  normalized  spherical  Bessel  equation  b(/^)  -  0  and 
=  0,  Substituting  (3.20)  into  (3.18),  reducing  by  means  of  the  differential 
equation  for  to  terms  containing  only  and  ,  and  equating  coefficients  of 
like  powers  of  /^p  and  on  both  sides  of  the  equation,  we  obtain  the  two  equations 

=  (3.21) 


ay(ay+l)-|.  ,  ZfeyfO  -  ^  rt  ' 

I  -  “5r*Py  . 


(3.22) 


These  are  two  simultaneoiis  equations  which  comprise  recnirrence  relations  for  the 
coefficients  Ay  and  By  in  (3.20).  If  we  choose  =  1,  then  the  first  few  coeffi¬ 


cients  are 


I 


mVp*-! 

■4a* 


26 


Thus  combining  (3.20)  and  (3.16)»  we  have  found  a  solution  of  (3«15)  of  the 

form 

u(z) -  -h  (3.23) 

where 

(3-24) 

®'*>  *  •  <3.25) 

In  order  to  conform  to  the  type  of  integral  representation  given  for  the 
spherical  case  in  See.  2,  we  choose  the  function  ^pto  be  the  normalized  spherical 
Hankel  function  h^  or  h^  .  Then  in  finding  the  normal  mode  solution  for  the 
spheroidal  problem  we  are  led  to  a  determination  of  the  zeros  of  the  function 

+  (3.26) 

We  can  reduce  this  problem  to  one  of  exactly  the  same  kind  as  solved  in 
Sec.  2.5.  From  (2.70),  we  can  replace  h^  (z)  by  a  suitable  sum  of  h2(^)  and  h^(i;) 
as  follows: 

h?^(a)  «  (3.27) 

From  this, 

(3.28) 

where  use  has  been  made  of  (2.62)  to  eliminate  hjgCi;). 

Introducing  (3.27)  and  (3.28)  into  (3.26),  we  obtain 

e  (3.29) 

where 

’S  (3.30) 

A  (S)  5  ai«W+  »  [4f5)t  (3.31) 

(3.29)  now  is  of  the  same  form  as  (2.70).  Consequently  the  procedure  by  which  the 

zeros  of  (2.70)  were  found  say  be  applied  directly  to  (3.29),  the  only  change 
required  being  the  replacement  of  and  jSC^)  by  a^^  (S)  and  jS,  ,  respectively. 
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4.  SI3MMARY 


In  this  report  ve  have  shovm  how  the  exact  earth-flattening  procedure, 
developed  in  [1]  for  an  Isotropic  apherlcally-stratlfied  a'haosphere,  may  be 
extended  to  the  ease  of  a  spherical  earth  and  alaiosidiere  enveloped  by  a  sharply 
bounded  ionosphere.  The  general  solution  of  the  problem  is  foxmtulated  as  an 
integral  representation,  which  nay  be  derived  either  a  ray-optical  series  or 
a  nonoaX  node  series.  In  the  latter  case,  the  normal  nodes  involve  the  nozmal- 
ized  spherical  Hjankel  function  and  its  derivative.  An  is^roved  method  of  obtain¬ 
ing  the  zeros  of  these  functions  is  derived  which  is  not  of  asymptotic  character. 

In  order  to  deal  with  problems  of  ncxi-spherical  stratification,  a  spheroidal 
gecmietry  is  investigated.  The  developments  for  the  spheroidal  case  are  pursued 
in  a  way  similar  to  that  for  the  spherical  geometry,  and  carried  out  in  detail  for 
the  oblate  spheroid.  Solutions  for  the  angular  function  are  found  in  the  form  of 
an  infinite  series  of  Bessel  functions  of  the  same  type  as  found  for  the  spherical 
case.  The  radial  function  is  expressed  as  a  sum  of  the  solution  of  the  normalized 
spherical  Bessel  equation  and  its  derivative,  the  coefficients  of  these  functions 
being  infinite  sexd.es  in  terms  of  powers  of  the  ratio  of  semi-focal  distance  to 
radius.  It  is  shown  that  the  zeros  of  the  radial  function  as  a  fmctlon  of  order, 
which  are  required  in  the  normal  mode  solution,  may  be  found  by  the  same  procedure 
that  was  developed  for  the  sphexdcal  case. 
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i^T  II 

VLF  ENH/LNCEMENTS  AND  HF  FADEODTS  DURING 
SnDDEN  IONOSPHERIC  DISTURB&MCES 

1.  INTRODUCTION 

One  of  the  speetactilar  phenomena  of  the  ionosphere  is  the  stdden 
ionospheric  disturbance  (SID),  which  drastically  affects  high-frequency  coHununi- 
cation  circuits o  This  phenomenon  was  first  reported  by  Mbgel  [1]*  and  later 
investigated  in  detail  by  Delling^  [2]«  Dellinger  sunmarized  the  vazlous  phenomena 
associated  with  the  SID  and  conclvtded  that  the  distuz'bance  must  be  caused  by  solar 
ultraviolet  radiation.  One  of  the  associated  phenomena  occurs  on  very  low  fre¬ 
quencies,  and  it  is  this  phenomenon  that  forms  the  subject  matter  of  the  present 
study. 

In  1936,  Bureau  and  Maire  [3]  reported  that  abrupt  short-^ave  fade-outs  (dmaoted 
SWF  hereafter)  usually  were  accompanied  by  simultaneous  sudden  increases  in  the 
strength  of  atmospherics  received  on  very  low  frequencies  (vlf).  They  reported 
that  atmospherics  from  all  directions  were  reinforced  simultaneously,  that  frequencies 
from  27  to  40  kc/s  showed  the  sudden  increase,  but  on  12  kc/a  the  effect  was  rarely 
observed.  Later,  Sudden  cmd  Ratcliff e  [4]  reported  that  measurements  at  Cambridge 
of  the  phase  of  the  abnormal  (horizontally-polarized)  component  of  the  downccmilng 
waves  fron  GBR  on  16  ke/s  showed  an  anomaly  at  times  of  h-f  fade-out.  They  concluded 
that  an  SID  "has  a  marked  effect  at  the  level  of  reflection  of  the  low-frequency 
waves  (70  km),  this  effect  being  most  evident  as  a  decrease  in  reflection  height 
of  the  waves".  They  did  not  observe  "say  clear  indication  of  a  change  in  reflected 
wave  amplitude  at  the  time  of  the  phase  anonalles"  (SPA).  Bureau  [5]  then  pointed 

^Nmbers  in  brackets  refer  to  the  corresponding  references  in  the  BibliograpI^ 
on  p.  52. 
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out  that  his  observations  on  the  sudden  enhanceaeiat  of  a too spherics  (SEA)  accompany¬ 
ing  SID  showed  that  such  increases  were  not  observed  belovr  aboui>  17  kc/s<. 

An  investigation  was  undertake!!  in  1938  to  determine  whether  SIDs  which  had 
been  shown  to  produce  SEA,  also  produced  sind-lar  enhancement  of  v-l-f  radio  signals, 
and,  if  so,  whether  any  quantitative  correlation  existed  between  the  v-l-f  end  h-f 
effects  of  the  SID.  The  experimental  phase  of  the  Investigation  was  completed  in 
1940,  and  a  preliminary  report  of  the  results  was  presented  in  3.947  [6],  but  has 
not  been  published. 

The  purpose  of  this  report  is  to  present  the  essential  results  obtained,  and 
to  discuss  the  implications  of  these  results  with  respect  to  ionospheric  layer 
structiire  and  the  modifications  produced  therein  by  the  SID  mechanism- 

2c  DESCRIPTIOM  OF  MEASUREMEHTS 

The  measurements  reported  here  were  macs  at  the  Riverhead  transcontinental 
receiving  station  of  RCA  CoBBiunicatlons,  Inc.  After  several  months'  observations 
of  the  signal  from  SAQ  (17.2  kc/s),  with  negative  rssxilts,  the  equipnent  was  set  up 
to  record  GLC  (31ol5  kc/s).  Some  of  the  subsequent  SWP'  were  accompanl-ed  by  sudden 
signal  enhanconents  (SSE)  of  GLC.  Consequently,,  ob^ejvratlons  were  continued,  extend¬ 
ing  over  the  period  31  October  1938  to  25  June  ’.940. 

For  comparison  of  the  v-l-f  SSE  with  SWF,  “he  signal  received  from  GLH  (13.53  Mc/s) 
was  selected,  since  this  signal  traversed  approximately  the  same  path,  and  continuous 
recording  of  this  signal  was  being  carried  out  fit  Riverhead  for  other  purposes.  The 
great  circle  path  length  was  about  5400  kaic  Bo'^b  the  GLC  and  the  GLH  equlmeats 
were  calibrated  at  least  once  each  day  Irr  means  o.f  standard  signal  generators = 

3.  RESULTS 

Sample  records  of  a  simultaneous  SWF  eaad  SSE  are  reproduced  in  FigSo  1 
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and  2f  respectively <>  These  records  are  rather  typical  of  the  data  obtained, 
although  the  magnitudes  of  the  signal  change  varied  rather  widely  from  one  event 
to  the  next.  In  general,  the  characteristic  behavior  was  a  rather  sharp  initial 
change,  followed  by  a  trough  (or  crest),  and  then  a  gradual  recovmryo  Invariably, 
the  recovery  was  more  rapid  for  the  h-f  signal. 

Fig.  3  shows  histograms  of  the  number  of  coincidences  between  SWF  of  GLH  and 
SSE  of  GLG  during  the  period  of  the  observations,  and  of  GLH  SWF  over  a  longer 
period.  Coincidences  were  observed  only  during  the  daylight  hours  when  the  h~f 
fades  were  more  nimerouso 

Flgo  4  shows  similar  histograms  of  the  number  of  GLH  fades  of  intensity  classi¬ 
fied  as  "moderate"  or  greater  and  (HiC  enhancements  which  occurred  during  the  same 
period  of  observation.  This  shows  a  high  degree  of  correlation,  so  that  the  proba¬ 
bility  of  a  v-l-f  enhancement  la  very  high  if  the  h-f  effect  is  pronouced. 

Fig.  5  represents  a  test  to  detexmine  whether  any  correlation  exists  between 
the  amplitude  ranges  of  the  v-l-f  and  h-f  signals  during  an  SID.  The  points  are 
plotted  with  the  Increase  in  GLG  signal  (in  decibels)  as  abscissa  and  the  correspond¬ 
ing  decrease  (in  decibels)  of  the  GLH  signal  as  ordinate.  Points  vdth  an  upward 
arrow  attached  correspond  to  cooqplete  fade-out  of  the  GLH  signal. 

Examination  of  Fig.  3  shows  that  in  no  case  was  the  GLC  Increase  as  great  as 
that  of  the  GLH  decrease,  and  that  no  evident  correlation  between  the  magnitudes 
of  the  two  effects  exists.  The  largest  GLC  increase  (14*1  db),  for  exai^ple,  was 
associated  with  only  a  moderate  fade  on  GLH.  Conversely,  the  deepest  fade  of  GLH 
(57.3  db)  was  accompanied  by  only  a  small  increase  (2.3  db)  on  GLG. 

4.  DISCUSSION 

In  the  years  since  the  observations  described  above  were  coiq>leted,  a 
considerable  body  of  information  has  accumulated  concerning  SID  effects,  solar 
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phenamena,  and  ionospheric  structure  <>  Observations  of  the  type  presented  above, 
however,  have  not  been  published  previously.  It  is  of  interest,  therefore,  to 
exanine  the  results  obtained  in  the  light  of  present-day  knowledge.  In  particular, 
it  appears  that  these  results  have  Important  liiq[>lications  on  the  type  of  solar 
event  which  causes  the  SID,  and  on  the  layer  structure  and  responsive  mechanisms 
in  the  upper  atmosphere. 

A  plausible  qrualitative  esiplanatio^  for  the  h>f  and  v-l-f  effects  was  advanced 
at  an  early  date:  The  h-f  waves  are  reflected  by  the  £-  and/or  F-layers;  absorption, 
however,  takes  place  mainly  in  the  intermediate  D-region.  V-l-f  waves,  on  the  other 
hand,  undergo  a  waveguide  type  of  pTOpagation  between  the  conducting  earth  and 
the  conducting  D-region,  the  attoiuatioo  depending  on  the  conductivity  of  the  guide 
"walls".  Since  an  enhancement  of  D->region  ionization  sbo\LLd  increase  the  "wall" 
conductivity,  this  will  reduce  the  attenuation  of  v-l-'f  waves,  but  will  give  rise 
to  tnc3reased  absorption  of  h-f  waves  passing  through  the  D-region. 

It  will  be  shown  below  that  the  above  qualitative  explanation  must  be  modified 
and  made  more  precise  in  order  to  fit  the  observations.  In  particular,  it  will 
appear  that  a  sharpening  of  the  lower  boundary  of  the  D-region  must  result  from  the 
flax^o  In  order  to  bring  this  out,  it  is  necessary  to  examine  the  absorption  and 
reflection  processes,  as  well  as  the  changes  in  ionospheric  layer  characteristics, 
which  take  place  as  a  resxilt  of  a  solar  flare. 

4.1  H-f  Effects 

Appleton  and  Plggott  [7j  have  made  a  comprehei:  slve  study  of  h-f  absorption 
at  vertical  incidence  dining  a  period  extending  over  a  sunspot  cycle.  They  found 
that  absorption  was  definitely  under  solar  control,  since  it  varied  in  a  regular 
manner  with  solar  zenith  angle.  They  showed  'Uiat  the  bulk  of  the  absorption  is  of 
the  non-deviative  type,  and  that  it  must  take  place  in  a  layer  below  the  reflecting 
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level  of  the  £<-region<i  Furthermorey  they  showed  that  the  absorbing  region  cannot 
be  mez^ly  the  lower  portion  of  the  S~regiony  but  must  be  an  Independent  ionized 
region t  which  they  identify  with  the  D-reglono 

s 

The  evidence  which  led  Appleton  and  Plggott  to  the  above  conclueiona  was 
obtained  tvotx  three  types  of  behavior  t 

(1)  The  diurnal  variations  of  absorption  for  two  different  frequencies^ 
one  of  which  is  reflected  by  the  £-layer  and  the  other  by  the  F-layer,  have  sub¬ 
stantially  the  sane  dependence  on  the  solar  zenith  angle  o 

(2)  For  a  frequoicy  whose  reflection  level  shifts  during  the  day  from 
the  F-layer  to  the  £-layer,  or  to  a  sporadic  £-layer,  the  absorption  is  the  same 
for  reflection  from  either  layer  (apart  from  the  period  when  the  frequency  is  in 
the  neighborhood  of  fS,  %d\en  additional  deviative  absorption  takes  place). 

(3)  The  variation  of  absorption  with  frequency  can  be  e^qplained  only  on 
the  assumption  that  the  same  medium  is  responsible  for  absorption  over  the  entire 


frequency  range. 

For  non-deviative  absorption  (i.e.,  in  a  region  where  the  refractive  index  is 
substantially  unity) ,  Appleton  [8]  gave  for  the  absorption  coefficient  tc  in  a  region 
of  ionization  density  N  and  collision  frequency  v,  under  conditions  where  the  quasi¬ 


longitudinal  approximation  holds » 


2it€*  NV 
me  <»>,.>* 


(1) 


where  is  the  magnitude  of  the  longitudinal  component  of  the  angular  gyro  fre¬ 


quency,  and  the  +  sign  is  for  the  ordinary  wave,  the  -  sign  for  the  extraordinary 


wave.  The  absorption  of  the  ordinary  wave  is  appreciably  less  than  that  of  the 
extraordinary  wave  when  u/u^  is  not  too  large,  so  that  it  is  the  ortiiziary  wave 
which  then  is  measured.  It  can  be  seen  that  the  dependence  of  ts  on  the  collision 
frequency  v  t^ds  to  a  proportionality  to  either  v  or  l/v,  depending  on  whether 
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Is  small  or  large  compared  with  (u  +  In  the  former  case,  the  integrated 

absorption  at  vertical  incidence  for  a  wave  which  penetrates  the  absorbing  region 
and  Is  reflected  (with  negligible  devlatlve  absorption)  at  a  higher  level  then  Is 
given  \jy  an  expression  of  the  form 

JWds  *•  Fi%),  (2) 

where  A  Is  a  constant  end  F(x)  is  a  function  of  the  solar  zenith  angle,  xt  which 
depends  on  the  rate  and  process  by  which  free  electrons  disappear  (e.go,  recombina¬ 
tion,  attachment).  Appleton  and  Plggott  showed  that  the  frequency  dependence  of 
the  total  absorption  (as  measxured  by  an  effective  reflecticm  coefficient)  is  in 
very  good  agremaent  with  (2).  This  is  shown  by  Fig.  6.  Thus  it  follows  that 
v^«(u  +  ux,)^  throughout  the  absorbing  region.  Appleton  and  Plggott  thus  placed 
an  upper  limit  for  v  of  2<>lo'Vsec  in  the  absorbing  D-region. 

Information  regaining  the  electron  production  and  ro&oval  processes  in  the 
absorbing  region  can  be  derived  from  a  study  of  the  dependence  of  absorption  on  the 
solar  zenith  angle  y.  In  partlctCLar,  the  theoretical  relation  shows  that  the 
effective  reflection  coefficient  p  depends  on  y  a  relation  of  the  form 

[logpl  <  (cosx)f^,  (3) 

where  n  depends  on  the  ionosphere  model.  For  a  Chapman  layer  (constant  scale  height 
and  recombination  coefficient),  n  =  1.5>  while  if  the  recombination  coefficient  is 
proportional  to  the  ambient  pressure,  n  =  1.0.  Nicolet  [9]  showed  that  a  region 
of  mounting  temperature  with  height  wo^ald  have  a  lower  value  of  n  than  one  of 
constant  temperature. 

The  experimental  values  cf  n  determined  by  Appletcn  and  Plggott  range  from 
about  0.4  to  1.1.  Taylor  [10]  foimd  values  from  0.7  to  1.30.  Furthermore,  Appleton 
and  Plggott  [7]  found  a  winter  anomaly,  the  absorption  in  winter  being  distinctly 
higher  than  for  the  same  zenith  angle  at  other*  seasons. 
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The  estperimental  values,  although  not  completely  understandable  on  the  basis 
of  present  theoretical  knowledge,  definitely  show  that  the  absorbing  layer  is  not 
of  the  Chapman  type  (for  which  n  -  1<>5},  snd  suggest  that  the  region  has  a  positive 
tsnperature  gradioit. 

The  above  studies  of  ionospheric  absorption  have  been  concerned  chiefly  with 
vertically  incident  waves  o  Since  the  path  Imogth  through  the  absorbing  region 
increases  as  the  secant  of  the  angle  of  incidence  on  the  absorbing  liy^er,  the  types 
of  variation  described  hold  substantially  for  an  oblique  path  of  constant  lengtho 
It  should  be  pointed  out  that  Appleton  and  Piggott’s  findings  relate  to  normal 
h-f  absorption,  and  that  the  height  region  wherein  the  additional  absorption  during 
SID  occurs  cannot  be  localised  fron  their  measurements o 
4o2  Y-l-f  Effects 

Although  the  main  features  of  h-f  absorption  are  fairly  well  underatood, 
this  is  not  the  case  for  v-l-f  waves »  The  requisite  theory  is  much  more  complicated, 
since  variations  in  the  properties  of  the  iiqx>rtant  regions  of  the  ionosphere  take 
place  in  a  distance  comparable  with  a  wavelength o  This  necessitates  full  wave 
theory,  which  is  made  complicated  by  ^he  anisotropy  of  the  medltaa.  An  analyticcd 
theory  has  been  worked  out  only  for  special  variations  of  electron  density  and  criti¬ 
cal  ft'equency  with  height,  and  thmi  only  for  the  case  of  a  vertical  magnetic  field 
or  of  vertical  propagation.  More  recently,  numerical  procedures  have  been  introduced 
to  handle  more  general  situations,  but  results  are  available  only  for  a  limited 
number  of  combinations  of  parameters. 

Our  present  knowledge  of  D-regicn  structure  has  been  promoted  fay  studies  of  the 
propagation  characteristics  of  v-l-f  waves.  These  characteristics  will  be  summarized 
here  in  order  to  provide  a  background  for  the  subsequmit  discussion  of  D-region 
mechanisms. 
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Although  some  measurements  of  layer  height  have  been  made  at  very  low 
fjrequencies  with  pulse  techniques  (BroMn  and  Watts  [11],  Hellivell  [12],  the 
Pennsylvania  State  Ihilversity  group  [13] )>  the  most  extensive  and  detailed 
studies  have  been  carried  out  on  e-v  tsransmissions,  principally  by  English  workers 
[14-22]*  These  measurements  have  been  made  at  various  distances  extending  out 
to  about  1000  km. 

The  principal  charactesristics  of  the  ionospherically  propagated  wave  (the 
so-called  "sky  wave")  are  its  phase,  amplitude,  and  polarization.  The  phase 
depends  on  the  length  of  the  transmission  path  and  the  height  of  reflection.  The 
apparent  height  of  reflection  is  deduced  from  observation  of  the  amplitude  pattern 
versus  distance  produced  by  interference  between  the  ground  and  sky  waves,  and 
also  by  measuring  the  phase  difference  between  gro\md  and  sky  waves  for  different 
frequencies.  Variations  in  reflection  height  with  time  can  be  deduced  from  measure¬ 
ments  of  the  phase  variation  of  the  sky  wave  at  a  given  receiving  point.  Vor  this 
purpose  the  sky  wave  is  isolated  from  the  gromd  wave  by  means  of  a  special  .uitenna 
arrangement.  Observations  of  the  change  in  phase  of  the  sky  wave  are  especially 
useful  in  testing  solar  control  of  the  reflecting  medium. 

Measiirements  at  a  frequency  of  16  kc/s,  for  example,  show  that  a  distinct 
change  in  the  character  of  the  sky  wave  takes  place  in  the  neighborhood  of  400  km, 
corresponding  to  an  angle  of  Incidence  on  the  ionosphere  of  about  65°.  Consequently 
it  will  be  convenient  to  discuss  the  short  and  long  distance  measurements  separate¬ 
ly,  and  then  the  modifications  observed  dtiring  SID. 

4*2.1  Short  Distance  Characteristics 

The  measurements  at  short  distances  may  be  s'cmmarized  as  follows; 
(a)  Reflection  Height 

Typical  results  of  the  phase  lag  of  the  sky  wave  relative  to 
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the  ground  wave  are  shown  in  Fig.  7..  The  height  of  reflection  shows  marked  solar 
control  during  the  day,  in  accordance  with  the  relation 

h  •  ho  I-  Att)  log  [Ch  CX)]  >  (4) 

where  h^  is  the  value  corresponding  to  x  ~  0»  and  Ch(x)  is  the  Chafanan  function, 
which  reduces  to  sec  x  X  Isas  than  about  85° «  An  average  value  of  hg  is 
73+2  ka.  If  reflection  took  place  fr«n  a  Chapman  layer,  the  slope  A(t)  of  the 
height  vso  log  iCh(x}]  curve  would  be  the  scale  height.  Fig.  8  shows  curves  of 
hQ  and  A(t}  at  16  kc/a  through  the  course  of  the  year.  The  apparent  heights  at 
noon  and  night  near  Cambridge,  England  are  shown  in  Fig.  9°  Values  of  A(t)  run 
around  6  km,  which  ia  a  reasonable  value  for  the  scale  height.  Consequently  this 
result  was  iised  for  some  time  to  infer  that  the  reflecting  layer  was  of  the  Chapman 
type.  On  30  kc/s,  however,  a  mean  value  is  5°5  ±  0.1  km,  on  43  kc/s,  4*8  +  0.1  km, 
and  at  70  kc/s  ai'ound  3  km,  with  greater  variability  at  the  higher  frequencies. 

This  variation  of  A(t),  howevmc,  is  not  explainable  on  the  basis  of  a  Chapman  layer. 
It  should  be  noted  that  the  descent  from  the  night-time  height  starts  at  a  time 
very  close  to  groimd  sunrise  at  the  midpath  point. 

(b)  Polarization 

For  short  distances  of  100-300  km,  nearly  all  observations 
show  that  the  sky  wave  on  all  frequencies  from  16-1^  kc/s  is  approximately  clrctilar- 
ly  polarized  with  a  left-handed  sense  of  rotation.  The  polarization  remains  the 
same  through  an  SID. 

(c)  Amplitude 

In  view  of  the  approximately  circular  polarization  of  the  sky 
wave,  the  components  P22  Pi2  tensor  reflection  coefficient  [see  Part  I, 

p.  11]  are  approxijsately  equal.  The  diurnal  variation  of  the  cooqxaient  p^,  called 
the  "conversion  coefficient",  is  shown  in  Fig.  10,  and  its  seasonal  variation  in 
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Pigo  11,  for  a  frequency  of  16  kc/s«  Fig.  12  shows  the  frequent  trend  of 
for  different  seasons. 

Figs.  13  and  14  show  the  diurnal  variation  of  pi2  on  16  and  70  ko/s,  respective¬ 
ly,  in  suner  and  winter.  It  is  seen  that  a  pre-sunrise  drop  and  post-sunset  rise 
in  sBqplltude  takes  place,  with  an  essentially  constant  level  during  the  day.  (The 
snail  ripples  in  the  winter  daytlne  curve  are  considered  as  probably  being  due  to 
a  two-hop  wave.)  The  drop  in  aiiq}litude  begins  at  a  solar  zenith  angle  of  close 
to  98®. 

It  is  evident  that  the  daily  anplltude  variation  is  distinctly  different  from 
the  dally  height  variation  at  short  distances. 

4o2.2  Long  Plntance  Charaoteristica 

The  characteristics  inferred  from  measurements  over  longer 
distances  will  be  aumarlzed  in  this  section.  These  principally  cover  distances 
of  about  400-9^0  km,  but  will  also  include  scmte  deductions  made  from  obaervaticms 
over  distances  of  several  thousands  of  kilometers.  These  have  been  derived  from 
fotir  sources;  (l)  16  kc/s  observations  at  540  km,  (2)  a  series  of  observations  over 
the  Decoa  navigation  chain  at  frequencies  from  70  to  about  130  kc/s,  and  distances 
up  to  950  km,  (3)  phase  variations  at  16  kc/s  and  lower  frequencies  in  connection 
with  basic  studies  of  navigation  systems,  and  (4)  observations  of  the  v-l-f  spectral 
characteristics  of  ataiospherics. 

(a)  Reflectl<«  Height 

The  reflection  heights  determined  from  the  ground  interference 
pattern  fit  in  with  a  reflection  height  of  70  +  2  km  at  midday,  with  no  apparent 
variation  of  height  %rith  frequency.  This  agrees  within  a  few  kilometers  with  the 
measurements  near  vertical  incidence. 

The  diurnal  variation  of  reflection  height  is  Illustrated  by  Fig.  15,  for  a 
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frequency  of  16  kc/s.  This  is  completely  different  from  the  diurnal  variation  at 
vertical  incidence  shovm  in  Fig«  7,  In  fact,  the  height  variation  is  very  much 
like  the  amplitude  variation  near  vertical  incidence  ahovn  in  Figo  10.,  Similar 
typea  of  variation  were  observed  at  higher  frequencies y  the  sunrise  drop  in  height 
being  substantially  complete  at  midpath  around  sunrise o  This  is  shown  in  Figo  16, 
for  which  it  was  assumed  that  the  nighttime  height  was  90  kmo 

Pierce  [23]  reported  a  nonaal  diur/ial  phase  variation  at  16  kc/s  of  200°  +  30° 
over  a  5200  km  path,  while  Casselman,  Heritage,  and  Tibbals  [24]  measured  a  diurnal 
change  of  about  350°  +  30°  at  12 «2  kc/s  over  a  4000  km  path., 

(b)  Polarisation 

Measurements  of  the  polarisation  of  the  sky  wave  shoved  this 
to  be  linear  at  about  43°  to  the  verticals  This  represents  a  change  from  the  short 
distance  meaaurements,  which  gave  the  polarisation  as  approximately  circular 

(c)  Amplitude 

The  reflection  coefficient  at  oblique  incidence  ia  found  to 
be  greater  than  at  vertical  incidenceo  For  16  kc/s,  Bain,  et  si  [19]  found  a 
value  of  0o27  at  aunmer  midday,  and  0^55  at  night,  compared  to  vertical  incidence 
values  of  OdS  and  0o50,  respectlvelyo  For  higher  frequencies,  Weekes  and  Stuart 
[21]  obtained  the  results  shovm  in  Fig„  17o  This  shows  an  increasing  reflection 
coefficiait  with  distance,  but  smaller  values  at  increasing  frequency c  Also,  an 
increase  of  about  2 si  takes  place  between  summer  and  winter o 

The  drop  in  amplitude  around  sxmidse  is  shown  in  Fig.,  IS,  This  is  similar  to 
the  behavior  of  the  reflection  height  shown  in  Fig^  15,  and  to  the  amplitude  be¬ 
havior  at  short  distances o  Again,  smaller  values  of  reflection  coefficient  are 
found  at  the  higher  frequencies » 

From  measurements  of  v-l-f  transmissions  on  available  frequencies  analyzed  by 


Eckersley  [25] »  combined  with  observations  of  the  spectrum  of  individual  atmos- 
pherics)  Ghapoan  and  Macarlo  [26]  deduced  the  attenuation  vs.  frequency  curve 
shown  in  Fig.  19*  This  shows  a  minimum  around  15  ko/s,  and  a  maximum  around  2  kc/s. 
4.2.3  gffijffftgti 

The  effects  of  SH)  associated  with  solar  flares  have  been 
observed  both  at  the  short  and  long  distances  used  to  obtain  the  results  discussed 
above.  In  general,  a  change  both  in  phaae  and  amplitude  of  the  sky  wave  is 
associated  with  an  SID.  The  change  in  phase  corresponds  to  a  decrease  in  reflection 
height.  This  change  in  phase  appears  to  be  a  very  sensitive  way  to  detect  flares 
[27]. 

Near  vertical  incidence,  the  decrease  in  reflection  hei^t  is  substantially 
the  sane  for  i^equeneies  in  the  range  16-135  kc/s.  This  is  illustrated  by  Fig.  20(a) . 
The  amplitude  near  vertical  incidence  suffers  a  decrease  during  an  SID,  the  change 
in  aiiq>lltude  being  greater  at  higher  frequencies,  as  shown  in  Fig.  20(b).  The 
relative  change  in  amplitude  is  roughly  proportlcmal  to  the  decrease  In  reflection 
height,  as  shown  by  Fig.  21  for  16  kc/s. 

The  above  characteristics,  observed  near  vertical  incidence,  undergo  a  drastic 
change  at  oblique  incidence  associated  with  the  longer  ranges  (>500  km).  The  phase 
change  associated  with  the  reductltm  in  height  of  reflection  decreases  with  increas¬ 
ing  frequency,  while  the  amplitude  increases  markedly.  The  amount  of  this  increase 
is  greater,  for  example,  at  70  kc/s  than  at  higher  frequencias.  Fig.  22  shows  an 
ezaiiq>le  of  the  relative  phase  and  amplitude  changes  observed  at  a  distance  of  about 
900  km  during  an  SID.  From  obaorvations  of  SEA,  it  appears  that  the  aaq>litude 
^jncrease  may  be  a  maximum  for  frequencies  around  30  kc/s. 

Piwce  [23]  showed  an  example  of  a  phase  advance  at  16  kc/s  over  a  5200  km  path 
during  an  SID.  This  SID,  of  importance  3,  accompanied  a  solar  flare  of  ixqx}rtance  2+, 
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A  phase  advance  of  100®  was  observedo  This  is  half  the  normal  diurnal  change,  or 
equivalent  to  a  reduction  in  height  of  reflection  of  about  9  kHo  No  amplitude 
change  was  observed,  however.  On  the  other  hand,  during  a  3-  SID,  aconapanying  a 
2  flare,  a  60  kc/s  signal  over  the  same  path  e3(perieDced  a  phase  advance  of  only 
70®,  cozreeponding  to  a  height  change  of  about  1.6  km,  while  the  amplitude  Increased 
considerably.  Fierce  suggested  that  the  primary  physical  phenomenon  produced 
the  SIO  might  be  a  steepening  of  the  ionization  gradient,  with  an  accompanying 
reduction  in  the  phase  lag  at  reflection. 

Gardner  [28]  and  Obayashi,  et  al  [29,30]  showed  that  an  SID  shifted  the 
.-squency  spectrum  of  atmospherics  upwards,  so  that  the  frequency  of  minimum  attenm- 
tlon  was  raised.  Also,  the  low-frequency  cutoff  of  the  ionospheric  waveguide  was 
raised,  corresponding  to  a  decrease  in  height  of  the  reflecting  region. 

To  simssarlze  the  SID  effects  observed  v->l-f  wave  propagation,  the  SID 
produces  a  reduction  in  reflection  height  and  a  change  in  amplitude  of  the  sky  wave. 
Near  vertical  incidence  the  reduction  In  reflection  height  appears  to  be  substan¬ 
tially  independent  of  frequent^,  while  the  amplitude  change  is  a  decrease.  The 
amount  of  this  decrease  Is  progressively  greater  at  higher  frequencies,  and  roiighly 
proportional  to  the  decrease  in  reflection  height.  At  100  kc/s  the  decrease  may 
be  by  a  factor  of  about  100.  At  oblique  incidence,  on  the  other  hand,  the  decrease 
In  reflection  height  Is  less  for  higher  frequencies,  while  the  aky  wave  amplitude 
increases  markedly.  This  Increase,  which  may  be  by  a  factor  of  5  or  more,  appears 
to  be  a  maximum  at  frequencies  around  30  kc/s,  and  becomes  less  for  higher  fre¬ 
quencies. 

4.2.4  Eclipse  Effects 

Observations  of  the  phase  of  the  sky  wave  on  16  kc/s  at  steep 
incidence  were  made  during  a  partial  solar  eclipse  ly  Bracewell  [31].  Although  the 
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greatest  ecllpaed  area  was  only  0,3  of  the  solar  disk,  a  definite  phase  ancoialy 
was  fovmd,  as  shown  in  Fig,  23,  The  foxn  oan  be  seen  to  agree  roughly  with  the 
shape  of  the  obscured  area  curve* 

Froa  this  result,  Braoewell  deduced  that  the  relaxation  tine  of  the  reflecting 
region  probably  did  not  exceed  6  ninutes.  Furthermore,  the  magnitude  of  the  phase 
change  —  about  35  degrees  —  represented  an  increase  in  height  of  reflection  of 
about  1  km,  while  for  a  Chapsam  layer  a  change  of  only  about  0,2  km  would  be 
expected , 

4o3  D-Layer  Production  and  Structure 

A  proper  interpretatioc  of  SID  effects  on  lonosi^eric  propagation 
ultimately  requires  a  knowledge  of  the  coaqpositlon  of  the  ionizing  agents,  and  of 
the  reactioas  which  lead  to  the  prevailing  ionization  densities.  In  this  Section, 
seme  of  the  pertinent  available  information  will  be  sunurized, 

40301  The  Two-Laver  Model 

In  order  to  explain  the  diurnal  phase  and  aaplitude  vazdations 
discussed  in  See.  4o2.1  and  4o2.2,  Bracewell  and  Bain  [32]  proposed  an  ionospheric 
model  contaixdng  a  two-layer  D-reglon.  The  height  of  the  upper  lay«:,  which  they 
denoted  by  Ikii  was  supposed  to  be  under  solar  control  in  accordance  with  the  formula 

h  =  72  +  5*5  log  sec  y  km,  (5) 

This  is  shown  by  the  upper  curve  in  Fig,  24 Below  this  layer,  a  layer  denoted  by 
was  postulated  to  exist,  with  height  variations  as  shown  in  the  lower  pert  of 
Fig,  24o  The  upper  layer  was  supposed  to  be  the  reflecting  layer  for  16  kc/s  waves 
at  steep  incidence,  while  the  lower  layer  was  consider^  to  be  respcmslble  for 
absorption  of  the  waves.  At  sufficiently  glancing  incidence,  however,  reflection 
would  take  place  at  the  lower  layer, 

Bracewell  and  Bain  baaed  their  two-layer  model  entirely  on  the  obsMrvations  of 


43 


16  kc/s  propagation  at  short  and  asdium  distances.  They  gave  no  suggestions  as 
to  the  mechanisms  by  which  these  two  layers  could  be  formed. 

4o3«2  Bracewell''a  EyhBustion  Region 

In  order  to  explain  the  observed  type  of  solar  flare  and  eclipse 
effects  on  the  D-region,  Braeewell  [31]  postulated  the  existence  of  a  so-<ialled 
"exhaustion  region" ,  in  which  the  ionizable  constituent  exists  in  a  small  concentra¬ 
tion.  With  respect  to  a  two-layer  D-region^  this  mechanism  was  supposed  to  take 
place  in  the  upper  region,  denoted  by  Da  in  Sec.  4°3ol. 

firaceuell  showed  that  an  exhaustion  region  would  explain  the  amount  of  change 
in  reflection  height  observed  during  a  partial  solar  eclipse,  whereas  a  much 
smaller  change  would  result  from  a  Chapmen  region.  Ke  also  showed  that  an  exhaustion 
region  would  produce  h-f  absorption  whose  variation  with  cosy  agreed  in  general  with 
experimental  observations. 

Braeewell  also  showed  that  the  characteristics  of  an  exhaustion  region  would 
explain  satisfactorily  the  observed  reductions  in  reflection  height  during  solar 
flares.  For  example,  a  reduction  of  15  km  in  height  would  require  an  increase  in 
intensity  of  the  incident  Ionizing  radiation  by  a  factor  of  15 «  However,  no  attempt 
was  made  to  deduce  the  accompanying  effect  on  the  amplitude  of  v-l-f  waves. 

4o3o3  Ionization  Mechanisms 

The  existence  of  several  separate  mechanisms  for  the  formation 
of  ionization  in  the  D-region  has  been  brought  out  in  the  last  decade  or'  so.  Brown 
and  Petrie  [32],  pursTiing  a  suggestion  attributed  to  Ratcliffe,  have  eva2.uated  the 
role  of  photodetachment  of  electrons  from  Og  ions.  This  ion,  formed  by  the  attach¬ 
ment  of  an  electrcm  to  a  neutral  oxygen  molecule,  starts  bidJLdlrig  up  in  concentration 
arcvuid  sunset,  resulting  in  the  disappearance  cf  the  normal  D-lsyer.  The  nighttime 
level  of  ionization  below  the  E-layer  is  malntaiaed  by  cosmic  rays,  which  nary  in 
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intensity  vflth  latitude «  Visible  light,  extending  down  into  the  infrared,  can 
supply  the  energy  required  to  break  up  the  attachment,  and  thus  liberate  free 
electrons 0  Since  visible  light  can  reach  the  altitudes  >50  km  appreciably  before 
ground  simrise,  electrons  released  by  the  photodetachment  process  build  up  D-layer 
ionization  appreciably  before  sunrise.  Brown  and  Petrie  [33]»  and  Ifoler  [34] 
showed  that  this  explained  satisfactorily  the  pre-sunrise  drop  in  amplitude  dis¬ 
cussed  in  Sec.  4*2.1.  Aiken  [35]  verified  the  fact  that  a  two-layer  D-region 
would  be  produced  at  sisirlse,  the  lower  layer  being  due  to  cosmic  rays,  and  the 
upper  layer  to  photoionization  of  nitric  oxide  by  Lyman-a  radiation.  Thus,  in  the 
two-layer  model  of  Bracewell  and  Bain  discussed  in  Sec.  4.3*1,  these  mechanisms 
would  account  for  the  layers  and  Da,  respectively. 

Nlcolet  and  Aikln  [36],  in  a  discussion  of  the  formation  of  the  D-reglon, 
pointed  out  the  following  mechanisms  of  ionization  which  are  possible  at  levels 
below  85  km: 

(1)  X-rays  of  X  <  10  A; 

(2)  lyman-a  radiation  (X  =  1215. 7A)} 

(3)  Ultraviolet  radiation,  X  >  1800  Aj 

(4)  Cosmic  rays; 

(5)  Photodetachment  by  visible  radiation. 

The  normal  S-layer,  which  is  ascribed  to  the  combined  affect  of  soft  X-rays  in 
the  range  30-100  A  and  ultraviolet  radiation  (Lyman-p)  is  penetrated  ly  cosmic  rays, 
ultraviolet  radiation  of  X  >  1^X3  A,  lyman-a  and  hard  X-rays  (X  <  10  A), 

Of  these,  cosmic  rays  and  hard  X-rays  are  capable  of  ionizing  all  almujspheric 
constituents.  In  addition,  Lyman-a,  due  to  a  narrow  window  in  Og  absorption  at  the 
Lyman-a  line,  can  penetrate  to  low  levels.  A  minor  constituent,  NO  (%  1  part  in  10^® 
was  propose;d  by  Nicolet  [37]  as  the  ionizable  constituent  responding  to  Lyman-a  to 


45 


accouat  for  the  daytime  D-layer.. 

In  view  of  the  presently-accepted  view  that  the  upper  part  of  the  D-region, 

Da,  is  due  to  photoionization  of  NO  l^man-Of  it  is  tempting  to  suppose  that  NO 
is  the  ionizable  constituent  iresponaible  for  the  eochaustion  region  postulated 
Bracevell*  The  concentration  of  NO  has  been  estimated  by  Nicolet  [38]  as  about 
10*”^®  of  the  total  concentiratioii  below  about  85  km,  or  about  10^  can-^  at  75  km 
[36].  In  order  to  give  ionization  densities  to  fit  the  v-l-f  obseirvations,  however, 
the  NO  concentration  wovtld  have  to  be  lower  than  this  by  about  two  orders  of  magni¬ 
tude,  or  about  10^  cm*"3  at  75  km. 

Although  Braceuell  believed  the  exhaustion  region  would  also  esqplain  solar 
flare  effects,  this  must  be  rejected  on  the  basis  of  later  evidence.  For  example, 
Friedman  and  collaborators  [39]  observed  no  large  Increases  in  lyman-a  during 
flares,  whereas  Bracewell  requires  a  factor  of  about  15°  In  a  recaat  report,  Chubb, 
et  al  [40]  stated  that  no  increase  in  lyman-a  occurred  during  a  1+  flare,  but  X-rays 
in  the  range  1-10  A  were  observed.  As  mentioned  earlier,  the  solar  flare  enhance¬ 
ment  of  ionization  has  been  shown  to  be  explainable  by  the  appearance  of  hard  X-rays 
in  the  wavelength  range  1-10  A,  which  ionize  all  atmospheric  constituents,  and  can 
penetrate  to  low  levels  because  of  the  low  absorption  coefficients  in  this  spectral 
region.  The  resulting  ionization  would  be  even  less  sharply  distributed  in  height 
than  a  Chapman  region. 

4.4  With  SID  Results 

The  two  features  of  the  experimental  results  shown  in  Fig.  5  which 
require  ejqplanation  are  the  followings 

(1)  The  lack  of  correlation  between  the  magnitudes  of  SWF  and  SSEj 

(2)  The  mechanian  which  produces  the  SSE. 

It  will  now  be  shown  that  the  first  is  explainable  on  the  basis  of  D-layer  structure 
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and  solar  flare  radiation,  but  that  an  adequate  explanation  of  the  second  is 
not  available  on  the  basis  of  present  knowledge. 

4.4ol  Afaiveace  of  Correlation  Between  Magnitudes  of  SWF  and  SSB 

The  absence  of  any  correlation  between  SWF  and  SSE  in  Fig.  5 
is  understandable  within  the  framework  of  the  two-layer  model  discussed  in 
Sec.  4°3<>1.  For  example,  if  the  flare  produces  a  burst  of  hard  X-rays  without 
any  enhancement  in  Lyman-a  radiation,  then  both  the  regions  of  the  Da  and 
layers  will  be  intensified.  The  relative  intensifications  of  these  two  regions 
will  depend  on  spectral  distribution  of  the  X-radlatlon.  There  is  no  reason 
to  believe,  at  present,  that  all  solar  flares  have  the  same  spectral  distributicn, 
so  that  the  relative  increases  can  be  expected  to  change  from  flare  to  flare.  The 
increase  in  h-f  absorption  leading  to  SWF  is  the  sum  of  the  increases  in  the  two 
regions,  \^le  the  v-l-f  SSE  would  respond  only  to  changes  in  the  lower  layer,  D^. 
Consequently,  this  would  result  in  the  absence  of  aay  clear-cut  statistical  correla¬ 
tion  between  the  v-l-f  and  h-f  effects  of  f]'.ares. 

4.4.2  Associated  With  SSE 

The  observations  reported  in  Sec.  3  show  that  SSE  on  vlf  is  one 
of  the  phenomena  accompanying  SID  {produced  by  solar  flares.  It  was  also  stated 
that  such  enhancements  can  be  understood  in  a  qualitative  way  as  due  to  reduced 
normalHsode  attenviation  as  a  result  of  Increased  conductivity  of  the  ionosphere, 
acting  as  the  upper  wall  of  a  waveguide.  It  will  now  be  sho^m  that  this  qualitative 
explanation  cannot  be  substantiated  on  the  basis  of  presently  accepted  ionization 
processes  and  present  theoretical  knowledge  concerning  v-l-f  propagation. 

For  the  ranges  involved  in  the  observaticns  reported  here,  the  nomalH&ode 
theory  of  propagation  is  more  advantageous  than  the  ray  theory,  since  only  one  mode 
is  effective.  A  nvnber  of  analytical  treatments  of  this  theory  have  appeared  [41-55], 
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but  none  treats  the  problem  in  a  sufficiently  general  way  to  definitiT© 

oonclueions  to  be  drawn  pertinent  to  the  present  observationss  “  ^iJaalytical  solutions 

have  been  obtained  only  for  speciial  distributions  of  ionizati  |®’^  density  and  colli »• 

’] 

sion  frequ«QQy  with  height,  and  f  £«*  special  dlrectloDS  of  th*||  j»a@aetlc 

field  (usually  taken  to  be  vertiosa).  Because  of  the  inabil  ,  ''^  produce  an 
analytical  solution  of  sufficient  generality,  efforts  have  directed  towards 
obtaining  numerical  solutions  [56“60j<.  This  approach  la  not  -^astrlcted  to  special 

height  distributions,  but  a  very  large  number  of  special  casj  ®  n®®d8  to  be  worked 

li 

out  in  order  to  produce  a  euff ici  imtly  extensive  catalog  frcsfij  deductions  of 

a  general  nature  can  be  dravna  in  yet,  only  a  rather  small  examples  has 

'j 

been  worked  out,  so  that  the  reavlto  from  which  one  must  drav'i  g®®®ral  inferwees 


are  rather  aoanty.  Nevertheleae,  these  tend  to  show  that,  otij^**  remaining 

unchanged,  the  attenuation  decreascc#  as  the  ionosphere  boundaiiy  sharper. 

Also,  for  a  constant  collision  fr'squonoy,  the  attenuation  deci  i*®®®*  ®® 
of  the  boundaxy  deoreaeee.  I 

One  of  the  idealizations  whl<3h  reduces  greatly  the  c<MQ5le3|^^y  ^  calcula¬ 
tions  is  that  of  a  sharply  bounded  homogeneous  ionosphero,  ca]j,|Culatioo8  using 
such  a  model  have  been  made,  amor,g  others,  by  Spies  and  Wait  ^der  the  further 

aa8unQ>tion  that  the  quaai-longiti-iSaaal  approximation  of  Booker  [j  ^^3  ®®y  ^  used. 

The  ionospheric  parameters  then  <?nter  the  analysis  in  an  8ffectj|r®  Wy 


given  by 


where  ujj,  v,  •r®  ■*^^®  pla«»f  collision,  and  longitudinal  gyrcl  »^S«lar  frequencies, 

respectively.  I 

Figo  25  (from  [53])  shows  the  attenuation  of  the  first  mode|i“  db/lOOO  km  as  a 
function  of  frequency  for  various  ionosphere  heights  for  a  value  “r  2°lo5„ 
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It  can  be  seen  fran  these  curves  that  a  reduction  of  height  from  75  to  60  km,  say, 
would  result  in  a  reduction  of  slightly  more  than  0o2  db/lOOO  km  for  a  frequenoy 
of  30  kc/so  For  a  5^\00  km  path,  the  total  reduction  in  attenuation  would  be  about 
lo2  db,  if  the  height  reduetion  ooourred  uniformly  over  the  ^ole  path.  This 
attenuation  decrease  is  the  result  of  a  decrease  in  the  grazing  angle  of  the  first 
mode  to  the  ionosphere o  Kwaver,  this  decrease  in  attenuation  is  based  on  a 
constant  effective  conductivity,  u^,  so  that  the  collision  frequen(7,  v,  is  assumed 
to  remain  unchanged . 

The  electron  density  distributions  in  the  D-region  shown  in  Fig.  26,  calculated 
by  Nioolet  and  Aiken  [36],  show  no  appreciable  change  in  shape  at  a  density  of 
about  10^  om*'^  between  a  quiet  sun  and  a  strong  flare.  Consequently,  a  solar  en¬ 
hancement  will  cause  a  given  ionizaticm  density  to  appear  at  a  lower  level,  but 
with  substantially  the  same  gradimit.  Hence  one  might  argue  that  it  is  reasonable 
to  suppose  that  a  decrease  in  attenuation  as  a  result  of  a  decrease  of  15  km  in 
reflection  height  of  the  same  order  as  that  calculated  for  the  sharply  bounded 
ionosphere  would  occur.  However,  in  virtue  of  the  approximately  e3qx>nential  ^crease 
in  critical  frequency  with  such  a  height  decrease,  the  value  of  u^.,  the  effective 
ionosphere  conductivity,  would  be  decreased.  On  the  basis  of  Kane’s  [62]  measurement 
of  collision  frequency,  a  15  km  height  decrease  would  bring  about  an  increase  in  v  of 
a  factor  of  10.  Assuming  a  value  of  ul  of  about  5'’10”  as  a  representative  value 
for  the  transatlant5.c  path  in  the  measurements  with  which  we  are  concerned,  then  the 
effective  conductivity  would  decrease  by  a  factor  of  about  2.3.  Thus  the  qualitative 
expectation  of  an  enhanced  ionospheric  conductivity  would  not  be  realized.  Instead, 
an  appreciable  decrease  in  effective  conductivity  would  resialt. 

The  above  conclusion,  it  must  be  emphasized,  is,  at  most,  semi-quantitative, 
since  it  is  based  on  the  behavior  of  an  idealized  sharply  bounded  ionosphere  having 
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"average"  properties  given  the  Nioolet  and  Aiken  results. 

In  order  to  obtain  an  increased  conductivity  at  the  lowered  heights  due  to 
the  onset  of  a  flarej  an  increased  gradient  at  these  lower  heights  appears  to  be 
required o  In  other  uordS)  in  addition  to  increasing  the  ionization  densities  at 
■n  levels  in  the  D*>regian»  it  appears  thiit  the  flare  mist  inorease  the  sharpness 
of  the  lower  boundazy*  This  would  result  in  a  decreased  penetration  of  the  waves 
reflected  therefron,  and  hence*  for  a  sufficiently  sharp  boundary*  could  outweigh 
the  effect  of  the  increased  collision  frequency  encountered  at  the  lowered  reflec¬ 
tion  height.  Again*  it  oust  be  enphasized  that  this  line  of  argument  is  only 
qualitative*  and  that  an  adequate  quantitative  theory  is  needed  before  a  fim  con¬ 
clusion  can  be  reached. 

If  we  grant*  for  the  tine  being*  that  an  increased  sharpness  of  the  lower 
boundazy  of  the  D-layer  is  required  to  explain  the  SSE  produced  by  the  flare*  then 
it  is  necessary  to  adduce  the  necfaanian  which  produces  this  effect.  As  mentioned 
above*  the  electron  density  distributions  calculated  by  Nicolet  and  Aiken*  which 
are  shown  in  Fig.  26*  show  no  appreciable  change  in  shape  at  the  electron  densities 
required. 

5.  CONCmsiONS 

Simultaneous  obsezvatloos  of  short-wave  fade-outs  (SWF)  of  a  13.5*^c/8 
signal  and  sudden  signal  enhancooents  (SSE)  of  a  31ol5>kc/s  signal  over  substantially 
the  sane  transatlantic  path  of  appz^3xinately  5400  kn  show  no  evident  correlation 
between  the  magnitudes  of  the  two  effects  of  the  SID.  This  absence  of  cozrelation 
is  understandable  on  the  basis  of  a  two-layer  D-region,  The  lower  layer  is  produced 
hy  cosmic  rays*  while  the  upper  layer  is  due  to  photoionization  of  nitric  oxide  bF 
Lyman-a  radiation.  Hard  X-rays  (in  the  range  1-10  A)  emitted  ly  a  solar  flare 
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penetrate  to  the  low  levels  of  the  D-reglon  and  ionize  all  constituents  (principally 
Og  and  Ng).,  The  relative  intensifioatiois  of  the  two  D-regicms  will  depend  on  the 
spectral  dlstrihution  of  the  X-radiation.  On  the  assutq)tion  that  the  spectral 
dlstribntian  varies  tram  flare  to  flare*  the  relative  Increases  also  can  be  expected 
to  vax7  fTon  flare  to  flare.  Since  the  increase  in  h-f  absorption  is  the  sus  of 
the  increases  in  the  two  regions*  while  the  v-l-f  tfibancesent  is  occasioned  only 
hy  the  changes  at  the  lower  level*  no  correlation  should  result  between  the  two 
effects. 

On  the  other  hand*  an  adequate  explanation  of  the  lasehaniem  of  the  v-l-f 
aihancesent  la  not  available  on  the  basis  of  present  knowledge.  Phase  aeasures^nts 
show  that  a  definite  decrease  in  height  of  the  lower  boundary  of  the  D-region  is 
caused  hf  the  flare.  This  reduced  height  causes  reflection  to  take  place  at  a 
level  of  higher  colllsi«xi  frequency*  which  should  result  in  a  decrease  in  the  effec¬ 
tive  conductivity  of  the  l^er  if  the  ionization  gradient  remains  the  saiee.  Conse¬ 
quently*  it  appears  that  an  increase  in  the  sharpness  of  the  lower  boundary  of  the 
D-region  is  required  during  the  onset  of  a  solar  flare.  The  mechanism  hy  which 
this  takes  place  needs  to  be  determined. 
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Fig,  1 _ Record  of  GLH  signal  showing  fade  due  to  SID  at  1510  U.T.,  March  20,  1939. 
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Fig.  2 — ^Record  of  GLC  signal  showing  enhancement  due  to  SID  at  1510  U.T.,  March  20,  1939 
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Tig.  24 —  The  diurnal  height  variations  of  tho  layers 
Da  and  Dp. 
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conductivity  Uj.  =  2x10^  for  various  heights. 
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